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Abstract 

We use the AdS/CFT correspondence to study propagation of sound waves in strongly 
coupled (2+1) dimensional conformal magnetic fluids. Our computation provides 
a nontrivial consistency check of the viscous magneto-hydrodynamics of Hartnoll- 
Kovtun-Miiller-Sachdev to leading order in the external field. Depending on the behav- 
ior of the magnetic field in the hydrodynamic limit, we show that it can lead to further 
attenuation of sound waves in the (2+1) dimensional conformal plasma, or reduce the 
speed of sound. We present both field theory and dual supergravity descriptions of 
these phenomena. While to the leading order in momenta the dispersion of the sound 
waves obtained from the dual supergravity description agrees with the one predicted 
from field theory, we find a discrepancy at higher order. This suggests that further 
corrections to HKMS magneto-hydrodynamics are necessary. 
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1 Introduction 

Gauge theory/string theory correspondence of Maldacena [1, 2, 3] proved to be a valu- 
able tool in study dynamical processes in strongly coupled four dimensional gauge 
theory plasma. One of the most impressive contribution of string theory to the non- 
equilibrium plasma phenomena was the construction of relativistic conformal hydro- 
dynamics [4, 5], with its emphasis on shortcomings of widely used second order hydro- 
dynamics of Miiller-Israel-Stewart (MIS) [6, 7]. 

It has been realized recently that string theory, and, in particular, dual holographic 
descriptions of strongly coupled (2+1) dimensional collective dynamics, might shed 
new light on longstanding problems in condensed matter such as the theory of quantum 
phase transitions [8], superfluidity [9, 10], and high-temperature superconductivity [11, 
12]. Our present paper is largely motivated by these ideas. Specifically, we would 
like to better understand conformal viscous hydrodynamics of strongly coupled (2+1) 
quantum field theories in external fields. Our starting point is the first order dissipative 
magneto-hydrodynamics proposed by Hartnoll et.al. (HKMS) in [8]. Much like four- 
dimensional MIS relativistic hydrodynamics, it is built on the idea of constructing 
an entropy current away from equilibrium plus constitutive relations for dissipative 
currents. The entropy current is constrained to have a positive divergence, an idea 
originally due to Landau and Lifshitz [13]. We already know that this framework 
misses some important aspects of relativistic hydrodynamics at second order [4], but 
it is a well-motivated approximation. In this paper we would like to subject HKMS 
magneto-hydrodynamics 1 to consistency tests by extracting transport properties from 
the dispersion relation of the sound waves in strongly coupled M2-brane plasma and 
comparing them with the transport coefficients obtained from analysis of the current- 
current correlation functions [15, 16]. 

Our second motivation is to use the external field as a 'dial' to control the (effective) 
speed of sound waves in strongly coupled plasma. Indeed, the speed of sound waves 
c s in conformal plasma and in the absence of external fields is determined by simple 
scale invariance, c s = l/\/2 in (2+1) conformal fluids, in units where the speed of light 
is unity. This result is universal in relativistic (2+1) dimensional conformal hydro- 
dynamics without external fields. On the other hand, interesting (2+1) dimensional 

1 We postpone detailed analysis of conformal viscous hydrodynamics in external fields to our com- 
panion paper [14]. 
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gapless systems, such as a single atomic layer of graphite [17], have a much smaller 
speed of propagating sound waves. Thus, if there is any hope of realistically modeling 
such systems in the context AdS/CFT duality, one needs to understand how to reduce 
the speed of sound in holographic relativistic plasma. In [18] it was argued that sound 
waves in 3+1 dimensions are coupled to the magnetic variables of the fluid since the 
magnetization oscillates with the fluid density. As a result, dispersion relation of the 
sound waves is affected by the external magnetic fields. We would like to understand 
here whether such an effect persists in strongly coupled (2+1) dimensional magnetic 
fluids which admit a dual holographic realization. 

In this paper, we explore viscous hydrodynamics of strongly coupled M2-brane 
plasma in external magnetic field. In the next section we discuss magneto-hydrody- 
namics from the gauge theory perspective, emphasizing the necessity to appropriately 
scale the external magnetic field in the hydrodynamic limit. We point out that de- 
coupling of the shear and sound modes in magnetic fluids requires vanishing of the 
equilibrium charge density (or the corresponding chemical potential). We derive dis- 
persion relation for the propagation of the sound waves, explicitly demonstrating its 
sensitivity to the background field, similar to what was observed in [18]. In particular, 
we show that we can reduce the speed of sound by appropriately turning on a mag- 
netic field. In section 3 we analyze magneto-hydrodynamics of dyonic black holes in 
AdS^ x S 7 supergravity backgrounds of M-theory. The latter realizes a holographic 
dual to strongly coupled M2-brane magnetized plasma. We study quasinormal modes 
of these dyonic black holes, and demonstrate, in parallel with the field-theoretic anal- 
ysis, that decoupling of the shear and the sound modes occur only for black holes with 
zero electric charge. We compute dispersion relation of sound quasinormal modes, and 
interpret the results within the framework of HKMS magneto-hydrodynamics. Finally, 
we conclude in section 4. 

2 Gauge Theory Magneto-Hydrodynamics 
2.1 Equations of Motion and Conformal Invariance 

In this paper, we are interested in hydrodynamic properties of the (2+l)-dimensional 
theory living on a large number of M2-branes. One can view this theory as the three- 
dimensional maximally supersymmetric gauge theory near the infrared fixed point. 
It also admits a holographic description as M-theory on the manifold AdS^ x S 7 [1]. 
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The state corresponding to the thermal equilibrium is described by the black brane in 
AdS± [19]. The most rigorous way to understand the field theory equations of motion is 
to use this holographic duality. For a review of the near-boundary holographic analysis 
of asymptotically AdS space-times see [20]. 

Let Qmn and A M , M,N = 0, . . . , 3 be the metric and the gauge field in AdS±. 2 
To construct the dual field theory on the boundary, we have to solve the equations of 
motion for Qmn and Am with appropriately defined boundary conditions gMN — ► gffl 
and A M -> \i,v = 0,1,2. According to the AdS/CFT dictionary [21, 19], the 
boundary correlation functions are encoded in the renormalized action S ren [g$, A^\. 
In constructing S ren [g$, A^] one can perform the integral over the AdS radial coordi- 
nate and it becomes a functional on the boundary. This procedure of constructing the 
renormalized boundary action is known by the name "holographic renormalization" . 
Then, the boundary stress-energy tensor and the current are defined as 

SSrenbUS, 4»>] = J ^xy/^g^T^Sg^ + (^f) . (2.1) 

The field theory equations of motion become the consequences of the symmetries of the 
action S ren . Since S ren is invariant under diffeomorphisms and gauge transformations 
it follows that we have the boundary conservation equations of the form 

V U T^ = F^J U , 

= 0, (2.2) 

where, to simplify notation, we have removed the brackets (). To derive eqs. (2.2) one 
has to assume that the diffeomorphisms and gauge transformations act non-trivially 
on the boundary. However, in the bulk there is a special type of diffeomorphism which 
does not transform the coordinates x M on the boundary and whose boundary effect is 
to Weyl rescale the metric g$ [22]. Such diffeomorphism can also act non-trivially 
on the boundary gauge potential [23]. In Appendix A we show that the induced 
transformation on the boundary gauge field is trivial for a massless bulk gauge field. 
Therefore, under such diffeomorphism, the boundary fields transform as 

Sg (P)i»> = -2a{x)g^ v , 8Af = , (2.3) 

2 For purposes of the paper, we will specialize to the case of the four-dimensional AdS space with 
three-dimensional boundary and restrict ourselves to the the metric and gauge bulk fields. Of course, 
this analysis can be made more general. See [20] and references therein for details. 
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where a(x) is the transformation parameter. See eq. (A. 4) in Appendix A. The renor- 
malized action S ren has to be invariant under all the bulk diffeomorphisms. The in- 
variance under these transformations implies that 3 

T/ = 0. (2.4) 

Therefore, conformal invariance is unbroken by the presence of the background gauge 
field. The unbroken conformal invariance can also be understood intuitively if one 
recalls that the fact that the gauge field goes to a finite x^ dependent piece A^} (x) on 
the boundary (see Appendix A), means that in field theory it represents a marginal 
deformation [19]. 

In this paper, we consider hydrodynamics in the presence of the net charge density p 
and the background magnetic field B. For simplicity, we will drop the label "(0)" from 
the boundary fields in what follows and denote the bulk and the boundary fields by 
the same letter. However, from the context it will be clear whether the corresponding 
field belongs to the bulk or to the boundary. 

2.2 The First Order Hydrodynamics 

In this subsection, we study hydrodynamic perturbations in the presence of the charge 
density p and the magnetic field B. For our purposes, we can take the boundary metric 
to be flat and the equations of motion become 

d v T^ = F^J U , 

<V = 0, (2.5) 

where F^ u is the field strength corresponding to the background electromagnetic field. 
One can do the standard decomposition of the stress tensor, 

= eu^u" + PA"" + IP" , (2.6) 

where 

A"" = if + u^u" , IPX = , (2.7) 

3 The right hand side of this equation has to be supplemented by the conformal anomaly which 
can also be computed holographically [24]. However, it is not relevant for our discussion since the 
conformal anomaly vanishes in odd dimensions. In the case of four-dimensional hydrodynamics it is 
relevant at fourth order as pointed out in [4] . 
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and e and P are the energy density and the pressure respectively. In writing eq. (2.6) 
we are working in the so-called "energy frame", and the four velocity of the fluid w M 
is defined by the eigenvalue equation T£u v = —eu^. It was argued in [8] that in the 
presence of the magnetic field the pressure P differs from the usual thermodynamic 
pressure p by the term —MB, where M is the magnetization. The dissipative term 
TV J,U is given by 

IF" = -na^ - (A^(d a u a ) , (2.8) 

where 

= A^A»P(d a u p + d p d a ) - A^(9 7 «t) > ( 2 - 9 ) 
and rj and ( are the shear and bulk viscosity. Note that IP 1 ' is, by definition, zero at 
local equilibrium. 

Now we consider the current J M . It is given by 

J" = fni" + v» , (2.10) 

where is the dissipative part satisfying u^v^ = 0. The expression for it can be 
obtained from the argument that the entropy production has to be positive. It was 
done in [8] and the result is 

^ = a Q A^(-d u p: + F va u a + ^d v T) . (2.11) 

In this expression, T is the temperature, \i is the chemical potential, F^ u is the back- 
ground field strength and oq is the DC conductivity coefficient. In the case under 
consideration 

F 0i = 0, i = 1,2, <,,P>. (2.12) 

We would like to study fluctuations around the equilibrium state in which 

^ = (1,0,0), T = const. , fi = const. . (2.13) 

As an independent set of variables we will choose the two components of the velocity 
Sui = Su x , 5u2 = Su y as well as 5T and 5/i. As usual, all perturbations are of the 
plane-wave form exp(—iuot + iqy). We find that the relevant fluctuations of are 

5T U = (e + P)5u l , 

V ; (2.14) 
ST* = -r)d y 5u x , 

5T yy = 6P _ {r] + C )d y 5u y = {^j T + (jp) ST-( V + ()d y 5u y . 
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Note that since P is different from the thermodynamic pressure, an d (§f)^ 

do not coincide with the charge density p and the entropy density s. The equalities, 
however, hold when we set B = 0. Similarly, we obtain the following fluctuations of 
the current 

5J X = p5u x + o-QBSuy , 

5J y = p5u y + a Q (-d y 5n + ^d y 5T - B5u x ) . (2.15) 

Substituting these expressions into equations of motion (2.5) and performing a Fourier 
transformation we get the following system of equations 



=u 



((|)/" + (^)/ r )- ,(£+P)fa - 



=w(e + P)5u y - q ( ( ^ 5p+ (j^pj 5T ) + iq\ v + ()5u y + ia Q B 2 5u y , 



+ iBpSu x , 

_ — n Rrr^( r) ii. — 

T 

I / r)n \ / r)d \ 

=00 



=u;(e + P)5u x — qBaq{5p — ^ST) — iBp5u y + iaQB 2 5u x + iq 2 r]5u x , 



(2.16) 

Our aim is to understand how sound and shear modes are modified in the presence 
of the charge density and the magnetic field. However, eqs. (2.16) are all coupled to 
each other and it does not seem to be meaningful to ask what happens to the sound 
and shear modes separately. On the other hand, there is a regime in which these 
equations simplify and decouple into the two independent pairs. First, we will consider 
hydrodynamics with no charge density and no chemical potential 

p = 0, p = 0. (2.17) 

In addition, motivated by M2-brane magneto-hydrodynamics, we will set 

(£X-°- (D/°< (2 - 18) 
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Conditions (2.18) are also satisfied at p = 0, fi = on the supergravity side which 
will be studied in the next section. Then, equations (2.16) get separated into the two 
decoupled pairs. The first pair reads 

( dP\ 

u{e + P)5u y - q ^— J 5T + tq 2 { V + ()5u y + tcT Q B 2 5u y = 0. (2.19) 

If we set B = 0, these two equations describe the sound mode with dispersion relation 

U j = ±c s q-tq 2 ^ 5 , (2.20) 

where the speed of sound is defined, as usual, as c 2 s = We will refer to this pair 
of equation as to the sound mode equations. The second decoupled pair of equations 
becomes 

uj(e + P)5u x — qBaqdfj, + iaqB 2 5u x + iq 2 r]5u x = , 



CO 



(jfj^J 5fi + qcr Q B5u x + iq 2 (T Q 5fi = 0. (2.21) 

In the absence of the magnetic field these two equations describe a shear perturbation 
5u x with dispersion relation 

oj = -^j^p , (2-22) 
and a diffusive mode 5ji with dispersion relation 

uj = -iq 2 , °\ . (2.23) 

\dn ) T 

We will refer to these equations as to the shear mode equations. 

In the presence of the magnetic field all these solutions (2.20), (2.22) and (2.23) 
disappear. From the sound mode equations we obtain one constant solution 

u = - % T^p + ° {q2) ' (2 ' 24) 

and a diffusive mode 

0J = - iq 2 cf-^ + O(q 4 ). (2.25) 

From the shear mode equations we also obtain a constant mode (2.24) and a subdiffusive 
mode 

B 2 (^) 



In this analysis it has been assumed that the magnetic field B is held fixed in the 
hydrodynamic limit. 

However, it is not clear whether the solutions (2.25) and (2.26) can be trusted. The 
first hint that the regime of constant B in the hydrodynamic limit might not be well- 
defined comes from inspecting the 5-dependent term in eq. (2.11). The conductivity 
coefficient oq is of order the free mean path I. Therefore, at B = each term in (2.11) 
is of order j where L is the scale over which the derivatives vary. However, the term 
with F ua at constant B is not of this order since this term does not contain derivatives. 
This means that this term is not small in the hydrodynamic limit | < 1. It also 
follows that the limit of small B does not commute with the hydrodynamic limit. In 
other words, B cannot be thought of as being a small perturbation and one can worry 
that the hydrodynamic analysis in this case is unstable under higher order corrections. 
Let us now present a more quantitative reason why the solutions given above might 
not be reliable. From equations (2.19) it follows that if uj ~ q 2 we obtain 

6T 1 (9 97) 

- . (2.27) 

ou y q 

Thus, assuming that the amplitude Su y is fixed and of order unity in the hydrodynamic 
limit, we find that the amplitude of 5T is infinitely large. This means that the terms 
which are naively of higher order because they are suppressed by higher power of u and 
q can, in fact, modify hydrodynamics at lower order because of the large amplitude. 

Hence, it is more natural to study magnetic fields which vanish in the hydrodynamic 
limit. That is, we consider B which scales as 

B = bq p , p > , (2.28) 

with b held fixed. An interesting observation is that choosing different values of p we 
can probe hydrodynamics in different regimes. In this paper, we will concentrate on 
the sound waves for p = 1 and p = 1/2. Below we will present the field theory results 
and in the next section, we will study the holographic dual description. 

Similar analysis can also be performed for the shear modes. We will not do it in 
the present paper. 

2.3 Sound Waves in Magnetic Field 

If the magnetic field vanishes in the hydrodynamic limit, to leading order in q we 
can consider the various transport coefficients and susceptibilities evaluated at B = 0. 
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This, of course, is consistent if we are interested in the first order hydrodynamics. If 
one wishes to go to the second (or higher) order one has to keep in mind that there 
will be corrections not only from the term which are of higher order in derivatives but 
also from the possible B = 6g p -dependence of the transport coefficients. Thus, the 
shear viscosity 77 and the conductivity oq can be taken to be equal to their values 
at B = 0. Since, as shown in subsection 2.1, our theory is conformally invariant the 
bulk viscosity ( vanishes. Furthermore, the pressure P in this case becomes the usual 
thermodynamic pressure p. Therefore, we have 



dP\ (dp_ 
df). = \df 



= s. (2.29) 



In addition, we have another well-known relation 

5e = c 2 s 5P, (2.30) 

where c 2 = 1/2. With these simplifications, we have the following sound mode equa- 
tions 

u)5e — q(e + P)5u y = , 

uj(e + P)5u y q5e + iq 2 i]5u y + ioQB 2 bu y = . (2-31) 

We would like to study solutions to these equations when B scales as bq and bq 1 ! 2 . 

2.3.1 The Regime B = bq 

Substituting B = bq into eqs. (2.31) we obtain 

u5e — q(e + P)5u y = , 

u(e + P)5u y - ]-q5e + iq 2 (r] + a Q b 2 )5u y = . (2.32) 

We see that the effect of the magnetic field is to shift the shear viscosity by the amount 
(jQb 2 . This means that the modified dispersion relation is 

1 iq 2 r] + o Q b 2 

^ 7?" ^ 7Tp~' (2 ' 33) 

We obtain a sound wave whose speed is 1/ \pl of the speed of light and with modified 
attenuation. 
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2.3.2 The Regime B = bq 1 ' 2 

Substituting B = bq 1/2 into eq. (2.31) we obtain 



u5e — q(e + P)5u y = , 



u(e + P)5u y — ig5e + iq 2 i]5u y + iaqb 2 q5u y = , (2.34) 

and the corresponding characteristic equation becomes 

uo 2 {e + P) + iu{ m 2 + a Q b 2 q) - l -q 2 {e + P) = . (2.35) 

Note that if we take uq and (e + P) to be equal to their values at B = we cannot 
trust the uq 2 term in this equation since it will be modified due to B 2 dependence of 
oq and (e + P). Then the solution becomes the following sound wave 









2(e + P) 2 








4(e + P) 2 



= - j 2(^9 ± ^ ( 1 " T7-^ + ) • (2-36) 

Note that the speed of this sound wave is different from l/y/2 and is decreased in 
presence of the magnetic field. 

Now let us compute corrections of order q 2 to the dispersion relation. Going back 
to eqs. (2.19) we see that to obtain all necessary terms to the given order we need to 
expand (e + P), |^,ctq to next-to-leading order in B 2 = b 2 q. Note that (rj + () 
already multiplies q 2 and, hence, can be taken in the limit of the zero magnetic field. 
Furthermore, from subsection 2.1 we know that the theory is conformally invariant. 
Hence, 

P=- r (2.37) 

where, as we have explained before, P = p — MB, with p being the thermodynamic 
pressure and M being the magnetization. Then it follows that eq. (2.35) still holds but 
we have to expand (e + P) and <jq to next-to-leading order in B 2 = b 2 q. Let us denote 

e + P = E + qb 2 E 1 , 

°q = (To + qb 2 (Ji ■ (2.38) 
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For simplicity, we will work to order 6 4 . Then we obtain the following solutions to (2.35) 



iq (a b ( E x . 4 



± 



9 A <rg& 4 r2 „ q . 



This finishes our field theory consideration. Now we are going to move to the 
supergravity side. Our aim will be to reproduce the solutions discussed in this section. 

3 Supergravity Magneto-Hydrodynamics 

3.1 Effective Action and Dyonic Black Hole Geometry 

The effective four-dimensional bulk action describing supergravity dual to M2-brane 
plasma in the external field is given by [15] 4 



1 < „4 



S A = — I dx 



--R+-h MN h -- 



(3.1) 



where 

1 = ^, (3.2) 

and gf is the bulk coupling constant. From eq. (3.1) we obtain the following equations 
of motion 

Rmn = ^F ML Fj^ — -gMNFhpF LP — 3qmn , 

1 (3.3) 

V M F MN = . 

According to the AdS / CFT dictionary, the equilibrium state of magneto- hydrodynamics 
is described by dyonic black hole geometry whose Hawking temperature is identified 
with the plasma temperature on the field theory side. A dyonic black hole in AdS^ 
with planar horizon is given by the following solution to eqs. (3.3) [15] 

ds\ = - Cl (rf dt 2 + c 2 (rf [dx 2 + dy 2 ] + c 3 (r) 2 dr 2 

= ^[-f(r)dt 2 + dx 2 + dy 2 ]+\^f- 
v t f(r) (3.4) 

F =ha 2 dx A dy + qadr A dt , 

f( r ) =1 + (h 2 + q 2 y -a + h 2 + q 2 y , 



^or simplicity, we set the radius of AdSn to unity. 
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where h, q, a are constants related to the field theory quantities as follows [15] 

R h 2 AnT a 2 B 2 

B = ha , fx — -qa, = 3 -, (3.5) 

a cr a 4 

where B is the external magnetic field of the equilibrium M2-brane plasma, p is its 
chemical potential, and T is the plasma temperature. 

Let us now review thermodynamics of dyonic black holes. It has been studied 
extensively in [15] to which we refer for additional details. In the grand canonical 
ensemble the thermodynamic potential is given by 

la 3 / u 2 B 2 \ 
n = -V 2P = V 2j2T (-l-£ + 3 _), (3.6) 

where V 2 is area of the (x, y)-plane and p is the thermodynamic pressure. Just like 
on the field theory side, as the independent quantities we consider a related to the 
temperature by eq. (3.5), the chemical potential p and the magnetic field B. In terms 
of these variables we can compute the energy e, the entropy s and the electric charge 
p per unit area. One obtains [15] 



J_ a 3 / B 
g 2 2 V a 2 a 



Just like on the field theory side, it coincides with the temporal component (Too) of 
the stress-energy tensor. Furthermore, the entropy density is given by 



Finally, the charge density is 



s = ^a 2 . (3.8) 

g 2 



p = \ ap. (3.9) 



In addition, we introduce the magnetization per unit area 

1 fdtt\ IB . . 

M = -—[ — ) = — --. 3.10 

In the presence of the magnetic field, the thermodynamic pressure p is different from 
(T xx ) by the term proportional to the magnetization. Just like on the field theory side, 
we introduce the pressure P as (T xx ) which equals to 

P = p-MB. (3.11) 
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It is easy to check that 

P = |. (3.12) 
Notice that the trace of the stress-energy tensor 

(T/) = (3.13) 

vanishes, implying unbroken scale invariance. 

In parallel to the field theory discussion in the previous section we study sound 
quasinormal modes of the magnetically charged black hole, i.e., we set 

q = // = 0. (3.14) 

Then from eqs. (3.6)-(3.12) we obtain 

/ de \ 1_ 2a 2 7r(3a 4 - B 2 ) ( de_ 

This is in exact agreement with our field theory conditions 
3.2 Fluctuations 

Now we study fluctuations in the background geometry 

9mn 9mn + h-MN , 
A M A M + au i 

where qmn and Am (dA = F) are the black brane background configuration (3.4), and 
{fiMN, «m} are the fluctuations. To proceed, it is convenient to choose the gauge 

To be consistent with the field theory side, we will take all the fluctuations to depend 
only on (t, y, r), i.e., we have a Z 2 parity symmetry along the x-axis. Strictly speaking, 
for this parity to be a symmetry the reflection of the x coordinate must be accompanied 
by the change of the 5-field orientation: 

Z 2 : x^-x k h^-h. (3.18) 
14 



0, 



(2.18). 



(3.15) 



(3.16) 



(3.19) 



At a linearized level, and for the vanishing chemical potential [i = 0, we find that 
the following sets of fluctuations decouple from each other 

Z 2 - even : {h a , h ty , h xx , h yy ; a x } , 
Z 2 - odd : {h tx , h xy ; a t , a y } . 

Notice that naively the gauge potential fluctuation a x is parity-odd, while the other two 
components are parity-even. This is misleading, as it will turn out that au oc h, and 
thus a x is parity-even, while {a t , a y } are parity-odd. The first set of fluctuations is a 
holographic dual to the sound waves in the M2-brane plasma in the external magnetic 
field, which is of interest here. The second set describes the shear and diffusive modes. 
Note that if the electric charge of the black hole q and, hence, the chemical potential 
is not zero, is a linear combination of the terms some of which are proportional to 
h and some are proportional to q. Therefore, in this case au does not have a definite 
sign under parity (3.18). As the result, the two sets of fluctuations in eq. (3.19) no 
longer decouple. To say it differently, if both h and q are non-zero, the electromagnetic 
background (3.4) is not an eigenstate of the spacial parity x — > —x. However, the 
decoupling of the sound and shear fluctuations requires that it be an eigenstate. This 
is in complete agreement with our field theory analysis. Eqs. (2.16) decouple into the 
two separate pairs of equations only if we set /i — 0, p — 0. 
Let us introduce 

\2 u „—i<-<Jt+iq y „ f„\2 

-iujt+iq y 

=c 2 (r) 2 h xx = e~ M ™ c 2 (r) 2 H x 



h u =ci(r) 2 htt = e" 

\2 



\2 



h ty =c 2 (ry h ty = e-*"™ c 2 {rY H ty , 



h,, 



(3.20) 



hyy =c 2 (r) 2 h yy = e-^y c 2 (r) 2 H yy , 



-ie 



-iojt+iq y ~ 



0, 



where {H u , H ty , H xx , H yy , a x } are functions of the radial coordinate only and Ci(r) and 
c 2 (r) are defined in eq. (3.4). Expanding at a linearized level eqs. (3.3) using eqs. (3.16) 
and eqs. (3.20) we find the following coupled system of ODE's 



o =K + K 



In 



qc 2 

c 3 



+ tt [H xx + H, 



yy} 



In 



C2 
Cl 



2c\ 



c 2 

Q — (H tt + H xx ) 

C 2 



2 2 
C 3 u2 4 (tt , it \ , a C 3 u„2. 



+ uj 2 (H xx + H yy ) + 2ujq Hty) - 3 -f h 2 a 4 (H xx + H yy ) + 6 -f ha 2 q a x , 



(3.21) 



=H ty + H'ty 



In 



C1C3 



c c 
+ -| uq H xx - 4^ ha 2 [ha 2 H ty + u a x ) , 

c 2 C 2 



(3.22) 



15 



= H "x + 2 H 'xx 



In 



qc 2 



-3 J 



+ 2 



In* 

Cl 



+ 7^2 ( u (H X x — H yy ) — q -^(Hu + H xx ) 



- 2uqH ty - -§ h 2 a* (H xx + H yy ) + 2 -f a x , 



=/J ra + 2 H 'yy 



c, -I I 

In** 

<5 



In* 

Cl 



(3.23) 

1 (? ( c 2 
zc 1 \ c 2 



+ 2wg^ ) - -f (# ra + Hyy) + 2 -| /ia 2 g a x , 



=a" + a' 



In* 


' 4 ~ I 


C 3. 





LO 



^i 2 ) + 77% ha ' z [ <i( H tt + h xx + 



2c| 



1 yyv 



+ 2o; 



"77 -Hfa 



(3.24) 
(3.25) 



The number of the second order equations, of course, coincides with the number of the 
independent fluctuations. Additionally, there are three first order constraints associ- 
ated with the (partially) fixed diffeomorphism invariance 



=U I [H XX + Hyy}' + 



In 



C2 
Cl 



Hxx + Hyy))+q[ H' + 2 



In 



C2 
Cl 



H 



ty 



In 



H tt )+^uH' ty + Aha 2 ^, 



=q ^[H tt - Hj - 

c 2 ( 

=[lncic 2 ]' [H xx + Hyy}' - [lnc 2 ]' H' tt + -§ (u 2 (H xx + H yy ) + 2cug if, 



{H tt - H xx ) ) - 2°-\ h 2 a 4 {H xx + Hy y ) + £\ ha 2 q a 

The constraints are just the Einstein's equations obtained by varying the action with 
respect to the pure gauge metric components h tr , h yr and h rr . We explicitly verified 
that eqs. (3.21)-(3.25) are consistent with the constraints (3.26)-(3.28). 

Now we introduce the fluctuations invariant under the residual diffeomorphisms 
and gauge transformations preserving the gauge (3.17). Since we have five second- 
order equations and three constraints there must be two gauge invariant fluctuations. 
We find them to be 

q 2 c[ci 



(3.26) 
(3.27) 

(3.28) 



X 



7 JR =4^ H ty + 2 H 



yy 



2H T 



U) 2 d 2 C2 



q 2 c 2 
+ 2—— H u , 



LU* C. 



Za =a x + — ha 2 (H xx - H yy ) 
Zq 



(3.29) 
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Then from eqs. (3.21)-(3.25) and (3.26)-(3.28) we obtain two decoupled (gauge invari- 
ant) equations of motion for Z H and Z A 

=A H Z H + B H Z' H + C H Z H + D H Z' A + E H Z A , (3.30) 

=A A Z" A + B A Z' A + C A Z A + D A Z' H + E A Z H . (3.31) 

The connection coefficients {Ah, • • • , -EU} can we computed from (3.21)-(3.25), (3.26)- 
(3.28) and (3.29) using explicit expressions for the q's, see (3.4). Since these coefficients 
are very cumbersome we will not present them in the paper. 6 In the next subsection, 
we will present the explicit form of the equations (3.30) and (3.31) in the limit of small 
uj and q. This will be sufficient for our purposes. 

3.3 Boundary Conditions, Hydrodynamic Limit and the Sound Wave Dis- 
persion Relation 

3.3.1 Boundary Conditions 

According to the general prescription [25, 26], in order to obtain the dispersion relation 
(poles in the retarded Green's functions) we have to impose the following boundary 
conditions on the gauge invariant fluctuations {Z H , Z A }. 

• {Zh, Z a } must have incoming wave boundary conditions near the horizon (as 
r-1); 

• {Zh, Z a } must be normalizable near the boundary (as r — > 0). 

The second condition is imposed because coefficients of non-normalizable solutions 
appear as poles in the retarded Green's functions. Thus, setting them to zero will 
produce the dispersion relation. See [26] for details. 

Since the solution of interest is an incoming wave at the horizon it has the following 
general structure 

Zh{t) = [f(r)f z H (r) , Z A (r) = [f(r)f z A (r) , (3.32) 

5 To achieve the decoupling one has to use the background equations of motion, i.e., the decoupling 
occurs only on-shell. 

6 The precise form of the equations (3.30) and (3.31) is available from the authors upon request. 
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where f(r) is given by eq. (3.4) and the functions zh{t) and z A (r) are non-singular at 
the horizon. Then denoting 

lim z H (r) -> zP + , lim z A (r) -> zf ^ , (3.33) 

i — > 1 _ 1 _ 

we find from eq. (3.30) and eq. (3.31) that as x = 1 — r — > + the following two 
equations must be satisfied 

= atv 2 (h 2 - 3)\h 2 q 2 - 3q 2 + 4m 2 ) 2 (4 / 3 1 2 + tv 2 ) z% ] x (1 + 0{x)) 
+ 8(/i 2 - 3) 4 (/i 2 q 2 - 3q 2 + 4to 2 ) 2 (-h 2 q 2 + 3q 2 + 16& - 4ra 2 )q^ ((3 - /i 2 ):r) Ml 
x (l + 0(x)) , 

(3.34) 

and 

= 2ahtv 2 {h 2 - 3)\h 2 q 2 - 3q 2 + 4ro 2 )(-/i 2 q 2 + 3q 2 + 8f3 1 )x z { ^ x (1 + O(x)) 
+ q(h 2 - 3)\h 2 q 2 - 3q 2 + 4tD 2 )> 2 + 4$) ((3 - h 2 )xf^ 1 zf x (1 + 0{x)) , 

(3.35) 

where we have defined ro = uj/(2itT) and q = q/(2irT). From eqs. (3.34) and (3.35), 
the existence of a nontrivial solution to (3.30) and (3.31) with incoming wave boundary 
conditions implies 

=atv 2 q((3 - h 2 )xf^{3 - h 2 f(q 2 h 2 + 4tt 2 - 3q 2 ) 4 
x j( m 2 + 4/3 2 )(m 2 + 4/3 2 ) + 0(:r) 
Eq. (3.36) suggests the following critical exponents: 

/3i = ft = -i|. (3.37) 

A detailed supergravity analysis of the solutions confirms that (3.37) is indeed the 
correct choice. 

Since (3.30) and (3.31) are second order ODE's, each of them has two independent 
solutions. Analyzing (3.30), (3.31) as r — * shows that normalizability of zh implies 
that 

^(r)=C(r 3 ) as r^O. (3.38) 

While both solutions of za are normalizable as r — > 0, requiring a fixed background 
magnetic field at the boundary implies that 

z A (r) = 0{r) as r^0. (3.39) 
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3.3.2 Hydrodynamic Limit with h oc q 

The aim of this subsection is to find the holographic solution dual to hydrodynamics 
with the magnetic field B ~ q. For this we have to study physical fluctuation equations 
(3.30) and (3.31), subject to the boundary conditions (3.37) and (3.38), (3.39) in the 
hydrodynamic approximation, ro — > 0, q — > with ^ and \ kept constant. 
To facilitate the hydrodynamic scaling we parametrize 

h=\q\H. (3.40) 

Furthermore, we parametrize the sound quasinormal mode dispersion relation as follows 

tv = c s q- t -(Tq 2 + T h h 2 ^+0 (q 3 , qh) , (3.41) 

with {c s , T, T h } kept fixed in the hydrodynamic scaling. Without loss of generality we 
can assume that T does not depend on H 2 , while T h = Y h (H 2 ). We will look for a 
solution as a power expansion in q and introduce 

z H = z H ,o + i q z HA + 0(q 2 ) , z A = h (z Afi + i q z AA + £>(q 2 )) . (3.42) 

Notice the 0(h) scaling of z A in the hydrodynamic limit. This is motivated by the field 
theory analysis. If p = 0, the x components of the current in (2.15) is proportional to 
the magnetic field. Since the equations for zh and z A are homogeneous we can rescale 
Zh{t) and z A (r) to make Zh{t) be equal to unity on the horizon. It is convenient to 
make a choice that the leading solution zh,o equals to unity and all the subleading 
contributions ZH,i(r), ■ ■ ■ vanish on the horizon. That is, we impose that 

lim Zh{t) = 1 , lim z H = 1 , lim zhi = . (3.43) 

r—*l— r— >1_ ' r—*l— ' 

To leading order in the hydrodynamic approximation, and subject to the appropri- 
ate boundary conditions, we find from (3.30) and (3.31) that 

=4 i0 - (^(9(4 C 2 + r 3 - 4) 2 + 64r 2 (r - l)(r 2 + r + l)(4c 2 - 3)H 2 ) 

x (r 2 + r + l)(r - l)rj ^9(4 - 4c 2 - r 3 )(4r 3 c 2 + 8c 2 - 8 + 4r 3 - 5r 6 ) 

+ 128r 2 (r - l) 2 (r 2 + r + l) 2 (4c 2 - 3)H 2 ^j z' Hfi + (^(9(4c 2 + r 3 - 4) 2 ( 3 ' 44 ) 

+ 64r 2 (r - l)(r 2 + r + l)(4c 2 - 3)H 2 )(r 2 + r + l)(r - 1) 
x 81 (4c 2 + r 3 -4)r 4 z H ,o , 
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and 

Solving (3.45) subject to the boundary conditions (3.43) implies that 

z H ,o = r 3 , c, = ±-j=. (3.46) 

Given (3.46), it is straightforward to verify that (3.44) is satisfied as well. Thus, we have 
solved the equations of motion to leading order in the hydrodynamic approximation. 

Now we move to next-to-leading order. Using (3.46), we find from (3.30) and (3.31) 
the following two equations 



=4,i 



128r 2 (r - l) 2 (r 2 + r + l) 2 H 2 - 9(r 3 - 2)(5r 6 - 6r 3 + 4) , 



rA 3 (r 2 + r + l)(r- 1) 



HA 



576i (r 3 - 2)H 2 r A . 81(r 3 -2)r 4 
+ z a, - A ~ 2 h I ,L T\ *h,i (3.47) 



+ 



aA 3 AU A 3 (r 2 + r + l)(r-l) 

9 v / 2(2// 2 (9r h - 8r 2 ) - 9 + 18r)(r 3 - 2)r 4 



A 3 (r 2 + r + l)(r-l) 
and 



+ 



„ r(64(r 3 -l)(r 3 + 2)F 2 -27r(r 3 -2) 2 ) , Qi r 2 a , 18i ra 

+ A 3 (r 3 - 1) Za >° + "at ^ ~ ~at~ ^ 

9i r v / 2a(8r /l i/ 2 (r 3 - 1) + 8r 3 r - r 6 - 8r) 



2A 3 (r 3 - 1) 



(3.48) 



where 

A 3 = 64r 2 (r - l)(r 2 + r + l)H 2 - 9(r 3 - 2) 2 . (3.49) 

It is difficult to solve eqs. (3.47) and (3.48) analytically for all values of H 2 . However, 
for our purposes it is enough to find a solution to order H 2 . Expanding Zh,i, %a,i and 
T h perturbatively in H 2 , 

zh,i = z H>lt0 + H 2 z Hxl + 0{H A ) , z Afi = z Afifi + H 2 z AA1 + 0{H A ) , 

(3.50) 

T h = r M + 0(H 2 ) , 

we find that the solution satisfying the boundary conditions stated in eqs. (3.38), (3.39) 
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and (3.43) is given by 



z h,i,o —0 , z h,i,i — , z a,o,i — , 

i V2a(9 ln(r 2 + r + 1) + 6^ arctan (^) - ttv^) 



£.4,0,0 — 



y^j v; (3.5i) 



72 

1-2' ift >°-g- 
Here the last line comes from requiring that vanishes on the boundary. 

5. 3. 3 Comparison with Field Theory 

To summarize the results, we have obtained the following dispersion relation 

to = ±-L q - 1 Q q 2 + ? tf 2 q 2 ) + O (q 3 , q 2 #) . (3.52) 

Let us compare it with the field theory counterpart (2.33). For this we will rewrite 
eq. (3.52) in notation of (2.33). First, we will recall that 

Second, we will express the variables (H, a) which are natural to use on the supergravity 
side in terms of (ft, T) . We get 

where 

a = -g- . (3.55) 

Note that the relation is non-linear because a also depends on B = bq. However, for 
our purposes it is sufficient to take it to leading order. Substituting eqs. (3.53)-(3.55) 
into eq. (3.52) we obtain (up to higher order terms) 

1 if q 2 4 ft 2 q 2 \ , ocfi , 
to = —=q-- 1-?-^ + t- ■ (3.56) 

Comparing it with eq. (2.33), first, we find that 

-^ = tV (3-57) 
e + P AixT V ; 
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This result agrees with earlier calculations of the shear viscosity rj in [27, 28]. Second, 
we obtain 



e + P 3ag " 

Recalling results for e and P from subsection 3.1 we conclude that 



(3.58) 



*Q = j 2 ■ (3.59) 

Thus, we reproduced the result for the conductivity coefficient (to leading order in the 
magnetic field) obtained earlier in [8, 16]. 

We see that we have a perfect agreement with our field theory analysis and with 
the earlier calculations of the transport coefficients. 

3.3.4 Hydrodynamic Limit with h oc ^/|q[ 

In this subsection, we will study the holographic solution dual to hydrodynamics with 
the magnetic field B ~ yfq. For this we will consider the physical fluctuation equa- 
tions (3.30) and (3.31), subject to the boundary conditions (3.37) and (3.38), (3.39) in 
the hydrodynamic approximation, ro — > 0, q — > with — and — 7= kept constant. 

The analysis will be similar to the one in subsection 3.3.2. To facilitate the hydro- 
dynamic scaling we parametrize 

h=^\o\H. (3.60) 
The sound quasinormal mode dispersion relation is parametrized as follows 

^A^ + M-^r^ + r, h^j + O (q 3 , q V, qh\ 

=r q + ir\ q 2 + (9(q 3 ), (3.61) 

r =AiJ 2 + c s , r 1 = -^r + r,^ 

with {A, c s , T, T/J kept fixed in the hydrodynamic scaling. We look for a solution as a 
series in q, 

z H = zh,o + i q z H ,i + q 2 z H ,* + C(q 3 ) , z A = h (z Afi + i q z A ,i + q 2 z A ,2 + C(q 3 )) , 

(3.62) 

and impose (without loss of generality) that 

lim zh{t) = 1 =>■ lim zh,o = 1 & lim zh,i = lim zh,2 = 0. (3.63) 

r— >1— ) — >1— ' ) — >1— ' r^l— ' 
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To leading order in the hydrodynamic approximation we obtain from (3.30): 

=4,o - ; 4,o , (3-64) 

which gives rise to the solution 

zh,o = r 3 . (3.65) 
Using this solution for zh,o, we find from (3.31) the following equation for za,o 

0-z" | r3 + 2 J | 9r ^ 2r o ~ ^ f3 6 6) 

U ~ Za >° + (r3 - l)r Za >° + 8(r3 - l)(4rg - 3)H 2 r ' l ^ bbj 

Solving (3.66) subject to the boundary conditions (3.39) gives 

ZA '° = 3 S(4q-°3) " 6 ^ arctan ^ " 9 ^ + *■ + !))■ ( 3 - 67 ) 

This finishes our consideration of the equations of motion to leading order. Note 
that we were not able to determine any of the transport coefficients in the dispersion 
relation (3.61). Moving on to next-to-leading order in the hydrodynamic approximation 
we find from (3.30) and (3.31) 

2 , 9ir 2 (r 3 + 41^-4) , i r 2 f 3 

-Z H>1 - - Z H>1 + 2r 2 a(r 3 _ fj *A,0 " 32(r 3 _ l)2 (4r 2 _ 3^2 l^ 4 " 81r 

+ A32i (r 3 - 2)iZVr + (162r 3 + 512# 4 r 6 - 972 - 512# 4 r 3 )r 2 
- 576i (r 3 - 2)H 2 r 2 T 3 + 648r 4 ) , 



(3.68) 



and 



0-z" I r3 + 2 z > | 3rg(4Tg-2-r3)a 



r ( r 3_i) ^ 16/7V 3 (r 3 -l)(4r2-3) 

2 A ,o(r);{r ,ri,i7};r 



+ 16r// 2 (r 3 -l)(4rg-3) Z - ff ' 1 + ' 



(3.69) 



where the source term J source is a linear functional of za,o also depending on the trans- 
port coefficients and the magnetic field , {T ,Ti, H}, as well as explicitly on r. 7 



7 The explicit expression for J SOU rce is extremely lengthy and cumbersome and we find it meaningless 
to put it in the paper. It is available from the authors upon request. 
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It is straightforward to analyze the asymptotic solution to (3.68) near the horizon, 
i.e., as x = 1 — r — > + . For generic values of T we find that zh,i has a simple pole 
and a logarithmic singularity near the horizon. Namely, 

zh,i =C {-3^ lux} + finite, (3.70) 

where 

C = 181^ + lQi H 2 T - 9 . (3.71) 
Regularity at the horizon implies that C = 0. It leads to the following solution 



T = c s + AH 2 = -i-H 2 ±-^Jl-—H*. (3.72) 

9 v^V 81 V ; 

Without loss of generality, we can assume that c s is independent of H 2 , while all the 
H 2 dependence resides in A. Hence, we obtain 



c s = ±-^, A± = -i - =f — X —i= 1- — # 4 I , (3.73) 



(3.74) 



where the signs ± are correlated in the above expressions. Given (3.73), a nonsingular 
function zh,i subject to the boundary conditions (3.38), (3.63) is uniquely specified: 

*■' = 27(32^-9,) X ((l + r + r 2 )(12V3a r cta„_ 

+ 181n(l + r + r 2 )) -2^- (36 + 2^)^ + r 2 ) - 27r 3 ^ . 
Now we move on to eq (3.69). Its general solution is of the form 

4i (I * (1 " e3) ^°" rCe(0 + <V,i) (3-75) 

where Ca,i,i is an arbitrary integration constant and 



Jsourceiy) — ^ s- 



z Afi (ry i {T ,T 1 ,H};r 



3rg(4rg-2-r 3 )q 
16# 2 r 3 (r 3 -l)(4r 2 -3) 



Since 



9i> 

+ 16riJ 2 (r 3 -l)(4r2-3)^' 1 



J source , aS 7" ► 0_)_ , 

r 

J source 5 X 1 7" ^ 0_i_ , 
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(3.76) 



(3.77) 



it is clear from (3.75) that for any value of I\ we can adjust the integration constant 
£4,1,1 to remove singularity of Za,i at the horizon, x — > + ; the second constant (ob- 
tained from integrating (3.75)) can be fixed to insure that za,i vanishes as r — * 0, see 
(3.39). Thus, we cannot determine Ti at this order. This is not surprising, given that 
to determine the leading order transport coefficient T we had to consider the 0(q) 
order in the hydrodynamic approximation of Zu- 
lu order to determine T\ we have to consider 0(q 2 ) in the hydrodynamic approxi- 
mation for zh- The equation of motion for zh,2 takes the following form 



n _,„ 2 9* rV + 4rg-4) 



za,o,z h ,i] {T ,T 1 ,H};r 



(3.78) 

where T source is a new source term. 8 Although we can not solve analytically for zh,2, it 
is straightforward to construct a power series solution, first, for za,i and then for Zh,2 
near the horizon, x = 1 — r — > + . Much like in (3.70), we find that z H>2 is non-singular 
at the horizon, provided 

1 8„ 4| 2^(32^-45) 
4 9 9V162 - 64# 4 

where ± sign correlates with the corresponding sign in (3.72). Perturbatively in H 2 , 

r x = - \ =f iM# 2 + 0(# 4 ) • ( 3 - 8 °) 
4 9 

5. 5. 5 Comparison with Field Theory 

To summarize computations performed above, we have obtained the following sound 
wave 



We would like to compare it with field theory. First, we will consider terms of order q 
and compare them with the field theory counterpart (2.36). Following the same logic 
as in subsection 3.3.3 it is straightforward to check that terms of order q are in total 
agreement with eq. (2.36) if, as before, 

' °~Q — ~~o ■ (3-82) 



e + P 4ttT' w g 



8 Its explicit expression is too lengthy and cumbersome to be put in the paper. It is available from 
the authors upon request. 
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Thus, we have obtained an agreement with field theory at leading order in q. Note 
that it holds to all orders in H . 

Let us now compare terms of order q 2 in eq. (3.81) with the corresponding field the- 
ory prediction (2.39). It is easy to check (it has already been done in subsection 3.3.3) 
that the results agree for vanishing magnetic field. However, we have a disagreement 
at order q 2 H 2 . The field theory result at this order is 

a 2 H 2 

tt, ~ T ^/2?—- ) (3.83) 
9 

where eqs. (3.82) have been taken into account. On the contrary, the supergravity 
result is 

to~±^l q 2 H 2 . (3.84) 
9 

We have been unable to identify the source of this disagreement in the framework of 
HKMS viscous magneto-hydrodynamics. It is natural to expect that introduction of 
additional transport coefficients would resolve this puzzle. We hope to discuss this 
issue in more details in [14]. 



4 Conclusion 

In this paper we discussed propagation of the sound waves in magnetic fluids in (2+1) 
dimensions. We used the strongly coupled M2-brane plasma and the general set- 
ting of the holographic gauge theory/string theory duality to test relativistic viscous 
magneto-hydrodynamics of HKMS [8]. We found that HKMS magneto-hydrodynamics, 
generalizing the arguments of Landau and Lifshitz [13] in constructing dissipative en- 
tropy currents, adequately describes propagation of sound modes to the order in the 
hydrodynamic limit first sensitive to the external field. There is, however, a disagree- 
ment in next order in the hydrodynamic approximation (still in the context of the 
first-order hydrodynamics). Such a disagreement suggests that additional transport 
coefficients, beyond those introduced in [8], might be needed to describe dissipation in 
(2+1) magnetic fluids even to first order in the local velocity gradients. 

Finally, it would be interesting to analyze shear modes in (2+1) dimensional strongly 
coupled magnetic plasma. While the general theorem [29] guarantees that the shear 
viscosity attains its universal value [30, 31, 32], a deeper understanding of the prop- 
agation of the shear modes might help in constructing a general theory of relativistic 
viscous magneto-hydrodynamics. 
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A Bulk Diffeomorphisms and Boundary Weyl Transforma- 
tions 

In this section, we show that the holographic stress-energy tensor is indeed traceless 9 . 
For the sake of generality, we will work in a bulk of an arbitrary dimension d + 1 . We 
write the metric using the Fefferman- Graham coordinates [33] as follows 

ds 2 = \ (dr 2 + g^(r, x)dx»dx v ) . (A.l) 

The boundary condition on the metric is 

\img^(r,x)=g<$ , (A.2) 

where g$ is the background metric which the gauge theory is defined. We do not need 
to worry about the exact nature of the sub-leading terms. 

The bulk diffeomorphisms that preserve the form of the metric (A.l) satisfy the 
Killing's equations 

V r & = 0, V (M £r)=0. (A.3) 

We are interested in bulk diffeomorphisms that generate a Weyl transformation in the 
boundary metric g$ . One can show that such solutions have the asymptotic behavior 

£ r = ^ + ... , ^ = - l -d,a{x) + ... (A.4) 

It is straightforward to show that the transformation of the boundary metric under the 
diffeomorphism (A.4) is given by 

8gV> = lim (-£&„,) = 2*(x)gV> • (A.5) 
9 Up to a possible conformal anomaly, which is known to be zero in odd dimensions. 
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In order to compute the transformation of the boundary gauge field, we need to 
know the asymptotic behavior of the classical solutions of the bulk Maxwell's equations 



V M F MN = . (A.6) 
In the gauge A r = 0, Maxwell's equations take the form 

A, + (2 - d)A, + \g aP uK + g, a g aP A p + -^Lg^dp (V=gg^g a5 F, 5 ) = o , 
^ v g 

(A.7) 

d, (y/^g^Aj) = , (A.8) 



where / = rd r f. One can then look for asymptotic solutions to these equations in the 
form A^ ~ r k A$ (x) + One finds that the exponent k must satisfy 

k(2-d + k) = . (A.9) 

Therefore, the leading term near the boundary is finite (k = 0) 

lim^(r,a;) = 4 )(x), (A.10) 

and we associate it with the boundary gauge field. 

For the sake of completeness, let us mention that in the case of a massive bulk 
gauge field, the RHS of the first Maxwell equation (A.7) equals m 2 A^. Therefore the 
asymptotic behavior of the field is modified to 



r->0 

where 



lim^(r,x) = r k Af(x) , (A. 11) 



k= 2 



1 r 



d - 2 - y/(d - 2) 2 + 4m 2 . (A.12) 



We still identify A ( °\ the background gauge field in the field theory. 

We are now ready to study the transformation properties of the bulk gauge field 
under the diffeomorphisms (A. 4). Our gauge is A r = 0. However, one needs to check 
that the diffeomorphisms (A. 4) respect such gauge choice. That is, in general one 
might need to make a compensating gauge transformation to ensure 5A r = 0. The 
leading change in the A r component is 

lim5 f ,4 r = lim(-£M r ) = r k+1 A^d^a(x) . (A. 13) 

r— >0 r—*0 " 
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We see that the compensating gauge transformation must be subleading at the bound- 
ary. More precisely, we take 5\A M = &m\ and from the requirement (6% + 5\)A r = 
we find that 

limA = --L- r k+2 A^d^a(x) . (A. 14) 

r^O k + 2 p 

Finally, we find that the total asymptotic transformation of the gauge field in the 
directions is 

lim(^ + 5 x )A fl = -ka{x)Af{x)r k , (A. 15) 

which translates to 

5Af = -ka{x)Af{x). (A. 16) 

Inserting the transformation rules (A. 5) and (A. 16) into the boundary action (2.1), 
we obtain the Ward identity 

(T£) = k(J»)Af , (A.17) 

where k is given in eq. (A. 12) for a massive gauge field. For a massless gauge field 
(k — 0), which is the case of interest in this paper, we obtain a traceless stress-energy 
tensor. 
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